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ABSTRACT 

The recent analysis of observations taken with the EIS instrument on Hinode suggests that well constrained 
measurements of the temperature distribution in solar active regions can finally be made. Such measurements 
are critical for constraining theories of coronal heating. Past analysis, however, has suffered from limited sam- 
ple sizes and large uncertainties at temperatures between 5 and 10 MK. Here we present a systematic study 
of the differential emission measure distribution in 15 active region cores. We focus on measurements in the 
"inter-moss" region, that is, the region between the loop footpoints, where the observations are easier to inter- 
pret. To reduce the uncertainties at the highest temperatures we present a new method for isolating the Fe XVIII 
emission in the AlA/SDO 94 A channel. The resulting differential emission measure distributions confirm our 
previous analysis showing that the temperature distribution in an active region core is often strongly peaked 
near 4 MK. We characterize the properties of the emission distribution as a function of the total unsigned mag- 
netic flux. We find that the amount of high temperature emission in the active region core is correlated with the 
total unsigned magnetic flux, while the emission at lower temperatures, in contrast, is inversely related. These 
results provide compelling evidence that high temperature active region emission is often close to equilibrium, 
although weaker active regions may be dominated by evolving million degree loops in the core. 

Subject headings: Sun: corona 



1. INTRODUCTION 

The distribution of temperatures in the solar atmosphere 
holds many important clues as to how the solar corona is 
heated. Coronal loops observed at temperatures near 1 MK, 
for example, often have very narrow temperature distributions 
(Aschwanden & Nightingale 2005; Tripathi et al. 2009; War- 
ren et al. 2008) and are evolving (e.g., Winebarger et al. 2003; 
Ugarte-Urra et al. 2009; Tripathi et al. 2010; Mulu-Moore 
et al. 2011), suggesting that these loops are far from equi- 
librium. Coronal emission at higher temperatures (~4MK) 
appears to behave differently. There is some evidence that 
the high temperature emission in the core of an active region 
is close to equilibrium (Winebarger et al. 201 1; Warren et al. 
201 1), suggesting that heating events must occur at high fre- 
quency to prevent loops from cooling. 

This difference in behavior between loops at different tem- 
peratures appears puzzling, but may be explained by recent 
work on wave models of coronal heating, van Ballegooijen 
et al. (2011) and Asgari-Targhi & van Ballegooijen (2012) 
have studied the dissipation of Alfven waves in the chromo- 
sphere and corona. The heating rate that they derive is highly 
localized at the loop footpoint and heating events occur at 
high frequency. This implies that short loops that are strongly 
heated will have high apex temperatures and, because heating 
events occur frequently, they will be close to equilibrium. For 
longer loops that are heated more weakly the apex tempera- 
ture will be lower. For such loops it is possible that no equilib- 
rium exists regardless of the frequency of heating events (e.g., 
Serio et al. 1981; Peter et al. 2012). This would give rise to 
evolving loops at lower temperatures. 

The observational evidence for equilibrium loops at high 
temperatures, however, is limited. Winebarger et al. (2011) 
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and Warren et al. (201 1) presented the emission measure anal- 
ysis for small areas in two active regions. In their analysis they 
find emission measure distributions that are strongly peaked, 
suggesting that loops are not evolving through a broad range 
of temperatures. Because of the limited sample size, it is un- 
clear how general these results are. Tripathi et al. (2011), 
for example, have found somewhat broader emission measure 
distributions for two other active regions. Viall & Klimchuk 
(201 1) have analyzed the temporal evolution of the emission 
in yet another active region and find evidence for evolving 
loops, even in the core. 

In this paper we present a more systematic survey of ac- 
tive region core emission. It is well known that the amount 
of high temperature emission scales with the total unsigned 
magnetic flux (e.g., Schrijver 1987) and we use this metric to 
parametrize the observed active regions. We have selected 15 
observations that span a wide range of magnetic flux values 
(10 21 -10 23 Mx). For each region we compute the differential 
emission measure (DEM) in the active region core using ob- 
servations from the EUV Imaging Spectrograph (EIS) and the 
Atmospheric Imaging Assembly (AIA). We focus on inten- 
sities measured in the "inter-moss" region, that is, the region 
between the loop footpoints where we are measuring the prop- 
erties near the loop apex. Measurements of the entire active 
region would potentially combine emission from both loop 
footpoints and cooling loops and would require full models of 
the active region to interpret. 

To better constrain the DEM at high temperatures we 
present a new method for isolating the Fe XVIII emission in 
the AIA 94 A channel. The results of this method have been 
calibrated against spectroscopic observations of Fe XVIII 
974.86 A (Teriaca et al. 2012). We find that for regions with 
appreciable magnetic flux the DEM in the active region core 
is strongly peaked near 4 MK consistent with our previous re- 
sults. For regions with weaker magnetic fields the amount 
of high temperature emission diminishes significantly and the 
DEM becomes broader, consistent with the analysis of Tri- 
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FIG. 1 . — AIA and HMI observations of solar active regions. The regions are presented in order of increasing total unsigned magnetic flux. Every image at 
a particular wavelength is displayed with the same scaling. The green boxes represent the regions selected to compute the emission measure distribution. The 
numbers in brackets are the fluxes given in Table 1. Data for regions 1—5 are shown here. 



pathi et al. (201 1) and Viall & Klimchuk (201 1). 

The observation of strongly peaked emission measure dis- 
tributions in active region cores casts strong doubts on the 
Parker nanoflare model of coronal heating (Parker 1988) as 
it is commonly interpreted. The very small spatial scales ex- 
pected for magnetic reconnection relative to the 1 " (725 km) 
resolution of current coronal instruments suggests that ob- 
served coronal loops should be composed of many unresolved 
threads that are various stages of heating and cooling. This 
implies that the observed temperature distributions should be 
broad (e.g., Cargill 1994; Klimchuk & Cargill 2001; Cargill 
& Klimchuk 2004). A survey of hydrodynamic simulations 
of coronal loops by Mulu-Moore et al. (2011) suggests that 



for nanoflare heating models the temperature distribution has 
a power-law index of about 2 or less (EM ~ T a ) Our analy- 
sis, in contrast, shows that the temperature distribution in the 
core of an active region is often strongly peaked, with a ~3- 
4. It is possible that magnetic reconnection in coronal loops 
behaves differently than has been assumed. It could, for ex- 
ample, occur more frequently or on larger spatial scales than 
previously imagined. At present, however, it appears that the 
wave heating scenario suggested by van Ballegooijen et al. 
(2011) is more easily reconciled with the available observa- 
tions. 

2. OBSERVATIONS 
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FIG. 2. — The same as Figure 1 but for regions 6-10 in Table 1. 



Our aim here is to investigate the temperature structure of 
solar active regions systematically. The observations of in- 
dividual emission lines from the EIS instrument provide de- 
tailed temperature diagnostics and introduces a strong con- 
straint on our analysis. EIS (Culhane et al. 2007; Korendyke 
et al. 2006) is a high spatial and spectral resolution imag- 
ing spectrograph. EIS observes two wavelength ranges, 171- 
212 A and 245-291 A, with a spectral resolution of about 
22 mA and a spatial resolution of about 1" per pixel. Solar 
images can be made by stepping the slit over a region of the 
Sun and taking an exposure at each position. 

Telemetry constraints generally limit the spatial and spec- 
tral coverage of an observation. These constraints necessitate 
the selection of a limited number of spectral windows in a 



raster and not all EIS observations include all of the poten- 
tially useful emission lines. We have designed several EIS ob- 
serving sequences that contain all of the emission lines needed 
to compute emission measure distributions. Of particular im- 
portance are the observation of emission lines from Ca XIV- 
Ca XVII, which constrain the analysis at temperatures above 
3 MK (Warren et al. 2008). These studies have been run fre- 
quently and we used summary images to manually review the 
available data and select a set of observations that appeared 
to span a wide range of solar conditions. Each EIS raster was 
processed in the usual way to remove the CCD pedestal and 
dark current, identify any defective pixels, and calibrate the 
data. Intensities were then determined for each emission line 
of interest at every spatial pixel using Gaussian fits to the line 
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FIG . 3 . — The same as Figure 2 but for regions 1 1-14 in Table 1 . 

profiles. 

For each observation we determined the NOAA coordinates 
for the active region of interest at the mid-point of the EIS 
raster. NOAA region numbers, times, and solar coordinates 
are given in Table 1 . There are EIS observations taken during 
the interval considered by Viall & Klimchuk (201 1), but they 
do not include several of the high temperature Ca lines and are 
not optimal for emission measure analysis. For completeness 
we have included an EIS raster from this time. The observa- 
tions analyzed by Tripathi et al. (201 1) pre-date the launch of 
SDO and are not included here. 

For each observation we obtained full-disk AIA images 
(Lemen et al. 2012) and full-disk Helioseismic and Magnetic 
Imager (HMI, Scherrer et al. 2012) magnetograms from the 
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Stanford JSOC data center. AIA is a set of multi-layer tele- 
scopes capable of imaging the Sun at high spatial resolution 
(0.6" pixels) and high cadence (typically 12 s). Images are 
available at 94, 131, 171, 193, 21 1, 304, and 335 A. AIA im- 
ages are also available at UV and visible wavelengths, but 
they are not used in this analysis. HMI also images the full 
Sun at high spatial resolution (0.5" pixels) and high cadence 
(typically 45 s). To simplify the data management we se- 
lected all of the data within 300 s of the raster mid-point (the 
dates and times given in Table 1), processed the images with 
the standard aia_prep software to provide a common plate 
scale, and averaged the images at each wavelength together. 
For the AIA EUV images we divide each image by the expo- 
sure time. We then extracted a 400" x 400" region centered on 
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TABLE 1 

A Survey of Solar Active Regions : 



Region 


NOAA 


Date 




Y ce „ 


Am 




hot 


QL 


EIS File 


1 


1082 


19-Jun-10 01:57:44 


-306.4 


439.3 


2.87(19) 


4.08(21) 


1.36(04) 


2.2 


eis.U.20100619.014433 


2 


1158 


12-Feb-ll 15:32:13 


-248.4 


-211.8 


3.04(19) 


4.22(21) 


1.88(04) 


2.7 


eis.11.201 10212.143019 


3 


1082 


21-Jun-10 01:46:37 


162.9 


405.2 


3.29(19) 


4.72(21) 


1.05(04) 


2.0 


eisJl.20100621.011541 


4 


1259 


25-Jul-ll 09:36:09 


224.7 


323.4 


3.98(19) 


5.68(21) 


7.81(03) 


2.0 


eisJl.201 10725.090513 


5 


1150 


31-Jan-ll 11:25:19 


-470.9 


-250.6 


5.17(19) 


7.50(21) 


1.47(05) 


2.2 


eisJl-20110131.102326 


6 


1147 


21-Jan-ll 14:10:50 


26.6 


476.5 


6.49(19) 


1.02(22) 


2.28(05) 


3.6 


eis.11.201 10121.133954 


7 


1243 


02-Jul-ll 03:38:08 


-299.0 


216.6 


6.79(19) 


1.09(22) 


6.18(04) 


2.9 


eis.11.201 10702.030712 


8 


1089 


23-Jul-10 15:03:07 


-363.4 


-453.6 


6.96(19) 


1.12(22) 


1.84(05) 


3.5 


eis.U.20100723.143210 


9 


1109 


29-Sep-10 23:51:36 


361.5 


261.5 


7.19(19) 


1.16(22) 


1.50(05) 


4.3 


eis.11.20100929.223226 


10 


1193 


19-Apr-ll 13:32:20 


36.3 


363.5 


7.94(19) 


1.35(22) 


2.50(05) 


3.3 


eis.11.201 10419.123027 


11 


1190 


11-Apr-ll 12:00:42 


-492.6 


281.0 


8.61(19) 


1.37(22) 


1.11(05) 


3.0 


eis.11.201 1041 1.105848 


12 


1271 


21-Aug-ll 12:25:42 


-50.8 


150.8 


9.59(19) 


1.57(22) 


1.47(05) 


3.6 


eis.11.201 1082 1.1 05251 


13 


1190 


15-Apr-ll 01:17:19 


218.1 


304.4 


1.04(20) 


1.78(22) 


5.22(05) 


3.7,3.3 


eis.11.201 10415.001526 


14 


1339 


08-Nov-ll 19:14:27 


88.1 


258.4 


1.41(20) 


2.60(22) 


6.08(05) 


4.8 


eis.11.20111108.181234 


15 


1339 


10-Nov-ll 11:33:19 


406.0 


266.8 


1.48(20) 


2.73(22) 


9.75(05) 


3.7 


eis.11.201 1 11 10.100028 



a X cen and Y cen are the NOAA active region coordinates differentially rotated to the mid-point of the EIS raster. Am is 

the total area occupied by pixels between 50 and 500 G in cm 2 . $m is the total unsigned magnetic flux in Mx. I^ot is 

the total AIA Fe XVIII intensity in the field of view in DN s _1 . The parameter a is the slope of the emission measure 
distribution between log T 6.0 and 6.6. 




Total Unsigned Flux (Mx) 

FIG. 4. — The total AIA Fe XVIII intensity (I hot in Table 1) as a function 
of the total unsigned magnetic flux (<E> m)- The solid line is a power-law fit to 
the data. The total amount of high temperature emission in an active region 
varies strongly with the total unsigned magnetic flux. 

the NOAA active region coordinates. For each magnetogram 
the line of sight magnetic field is corrected to the radial value 
by dividing by the cosine of the helospheric angle. Represen- 
tative images are shown in Figures 1-3. For each observation 
we coaligned the EIS Fe XII 195.1 19 A raster with the AIA 
193 A image. 

For each HMI magnetogram we compute the total unsigned 
flux (<&m) for radial magnetic field strengths between 50 and 
500 G. The lower bound excludes the quiet Sun and the upper 
bound excludes sunspots. These limits were used by War- 
ren & Winebarger (2006) to study the relationship between 
the total unsigned flux and the total soft X-ray intensity. As 
in previous studies, they found a power-law relationship be- 
tween the total intensity and the magnetic flux, I sxr ~ $> b M , 
with b « 1.6. The values for the total unsigned magnetic flux 
we find here are similar to those from our earlier study, which 
used magnetogram data from the Michelson Doppler Imager 
(MDI) instrument on SoHO (Scherrer et al. 1995). Note that 
the absolute magnetic fluxes measured with HMI need to be 
reduced by a factor of 1 .4 to agree with those measured with 
MDI (Liu et al. 2012). 

The 94 A channel on AIA contains the Fe XVIII 93.92 A 
line, which is one of the most intense Fe XVIII transitions 
(e.g., Desai et al. 2005). Since this line is formed at a high 
temperature (7.1 MK) we expect that its integrated intensity 
will have a dependence on $m similar to that of the soft X- 



ray emission. Unfortunately, as can be seen in Figures 1- 
3 this wavelength range also contains emission lines formed 
at lower temperatures. The atomic data for this wavelength 
range is incomplete (O'Dwyer et al. 2010; Testa & Reale 
2012), further complicating the quantitative use of this chan- 
nel. As we describe in detail in Appendix A, it is possible to 
use a combination of AIA 171 and 193 A images to estimate 
the amount of contaminating emission in the 94 A channel 
empirically. Subtracting the estimated warm emission from 
the observed 94 A image isolates the Fe XVIII 93.92 A contri- 
bution. We have applied our algorithm to each active region 
observation of interest and the results are shown in the final 
columns of Figures 1-3. Comparisons with spectroscopic ob- 
servations of Fe XVIII 974.86 A are given in Teriaca et al. 
(2012). 

For each AIA Fe XVIII image we have computed the to- 
tal intensity in the active region above 2DN s _1 , which we 
consider to be the noise level introduced by the subtraction 
method, and list these values in Table 1 . In Figure 4 we show 
a plot of the total intensity as a function of total unsigned mag- 
netic flux. We also show a power law fit to the data and obtain 
a power-law index of 2.3, somewhat higher than our previous 
result using soft X-ray images. This exercise confirms that 
while our set of active region observations is not large, it does 
sample a wide range of solar conditions. The range of total 
Fe XVIII intensity varies by almost 2 orders of magnitude, 
froml0 4 tol0 6 DNs- 1 . 

The next step in this analysis is to manually select a small 
"inter-moss" region for each active region. These sub-regions 
were chosen if they were bright in AIA Fe XVIII but did 
not contain significant footpoint (moss) in AIA 171 A. The 
term moss refers to the footpoint emission of high tempera- 
ture loops which appears bright in emission lines formed near 
1 MK (see Berger et al. 1999 and references therein). We also 
attempted to avoid 171 A loop emission in the core of the ac- 
tive region, but for some observations this was not possible. 
These selections are indicated by the boxes display in Fig- 
ures 1-3. Note that the inter-moss region considered here for 
the 2010 July 23 active region is slightly different than that 
analyzed in Warren et al. (2011). The highest AIA Fe XVIII 
intensities are seen in the 2011 April 15 active region as a 
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FIG. 5 . — Example EIS active region observations. Selected rasters and line profiles are shown for the 15 April 201 1 active region. The profiles from the intense 
"bar" of emission are displayed in the lower panels. The profiles from the larger box are displayed in the top panels. The red curves indicate the Ca lines of 
interest. The blue curves indicate emission lines formed at lower temperatures. Below each profile the difference between the fit and the observed line profile are 
displayed. The vertical line indicates the l-cr error for the intensity in the spectral bin. 



bright "bar" of emission. For this active region we select two 
fields of view, one on the bright bar and another where the 
intensities are weaker, but more similar to the intensities ob- 
served in the other active regions. 
For each inter-moss region we extracted all of the relevant 



EIS data from each spectral window and averaged them to- 
gether to create high signal-to-noise line profiles. In comput- 
ing these averaged profiles missing data are not included. We 
then fit the line profiles with single Gaussians. The Ca XVII 
192.858A line is blended with Fe XI 192.813 A and a com- 
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FIG. 6. — The "inter-moss" emission measure distribution for 15 active regions. The field of view used to derived the intensities is indicated in each panel in 
Figures 1-3. For each measurement the emission measure distribution derived from MCMC is shown (thick red line) as well as the results from 250 Monte Carlo 
runs (black lines). The emission measure loci curves are color coded by ion. The slope of the distribution from 6.0 to 6.6 is indicated by the yellow line. The 
power-law index is also indicated on each plot. Regions 1-6 are shown here. 



plex of O V lines. We use the method outlined by Ko et al. 
(2009) to disentangle this blend. To insure consistency be- 
tween the fits to the Ca lines we use the widths measured for 
the Ca XIV 193.874A and Ca XV 200.972A lines to con- 
strain the fits to the other Ca lines. The width of the Ca XVI 
2008.604 A is set equal to that of Ca XV 200.972 A. The width 
of Ca XVII 192.858 A is limited to be within 0.05 mA of the 
width of Ca XIV 193. 874 A. Example rasters, line profiles, 
and fits are shown in Figure 5. 
The final EIS line list for each active region is generally the 



same used in Warren et al. (2011), except that we now add 
intensities for Ar XIV 194.396 A, a high first ionization po- 
tential (FIP) element that is useful for measuring the compo- 
sition. As before, we also include S X 264.233 A and S XIII 
256.686 A, but Ar is formed at a somewhat higher tempera- 
ture and has a higher FIP. Del Zanna (2012) have considered 
the relative intensities of some of these high FIP lines in a 
diffuse off-limb active region spectrum and suggested poten- 
tial problems with blends. Our intensities for these lines are 
approximately 50 larger and we are able to obtain consistent 
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FIG. 7. — The same as Figure 6 but for regions 7-12. 

results for these lines. The Ca XlV-Ca XVI lines are not avail- 
able for the 2010 June 19 active region, which is included here 
because it was studied by Viall & Klimchuk (201 1). 

The intensity that we observe with EIS is related to line 
emissivity and the emission measure distribution by the usual 
expression 

h = ^J zx(ne,T e )£(T e )dT e , (1) 

where e\(n e ,T e ) is the emissivity computed with the CHI- 
ANTI atomic database version 7 assuming coronal abun- 
dances (Feldman et al. 1992) and the CHIANTI ionization 



fractions (Dere et al. 2009). The function £(T e ) = n\ ds/dT e 
is the differential emission measure distribution and the chal- 
lenge we face is to infer this distribution from the observed 
intensities. It is also useful to consider the emission measure 
loci computed from 

£,loci{T e ) = j T^TTj (2) 

e\(n e ,T e ) 

which indicates the temperature range where the various lines 
are sensitive. Note that to aid in the comparisons with the em 
loci we will always plot the DEM multiplied by the tempera- 
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Fig. 8.— 
ture bin, 

aTe)dT e , (3) 

and we refer to this as the emission measure distribution 
(EM). 

We also wish to use AIA Fe XVIII intensities derived from 
our subtraction method to further constrain the emission mea- 
sure calculations. The hottest strong emission line observed 
with EIS during non-flaring conditions is Ca XVII 192.858 A, 
which is formed at about 5 MK. As mentioned previously, this 
line is blended, which adds considerable uncertainty to the 
intensity. To utilize the subtracted AIA Fe XVIII intensities 
we have computed a new response for this channel that only 
contains contributions from Fe XVIII and continuum. The 
response distributed with the official software contains con- 
tributions from several of the known emission lines formed at 
lower temperatures. 

To compute the differential emission measure we use the 
Monte Carlo Markov Chain (MCMC) emission measure al- 
gorithm (Kashyap & Drake 1998, 2000) distributed with the 
PINTofALE spectral analysis package. This algorithm has 
the advantage of not assuming a functional form for the dif- 
ferential emission measure. The MCMC algorithm also pro- 
vides for estimates of the error in the EM by calculating the 
emission measure using perturbed values for the intensities. 
The algorithm assumes the uncertainties in the intensities are 
uncorrelated so that systematic errors in the calibration, which 
could depend on the wavelength, or in the atomic data, which 
could vary by ion, are not accounted for. 
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For each "inter-moss" field of view we have run the MCMC 
algorithm to compute the differential emission measure. Ad- 
ditionally 250 Monte Carlo runs have also been performed 
for each field of view. The resulting temperature distributions 
are shown in Figures 6-8. The agreement between the ob- 
served intensities and those computed from the EM is gen- 
erally good, with most differences at the ±25% level. For 
some of the weakest active regions the intensities of the hot 
Ca lines become difficult to determine and the differences be- 
tween the observed and computed intensities are as much as 
50%. Inspection of quiet regions suggests that the Ca lines are 
all weakly blended. Quiet sun intensities are typically about 
5% of the intensity in the core of a very bright active region. 
For the weaker regions the impacts of the blends are more 
significant. We have not attempted to correct for these blends 
and so the observed intensities represent upper bounds. 

Inspection of the emission measure distributions indicates 
that many are strongly peaked near 4 MK (log T e = 6.6), sim- 
ilar to the result from our previous analysis (Winebarger et al. 
2011; Warren et al. 2011). To quantify the steepness of the 
emission measure we fit a power law of the form EM ~ T a 
to each distribution. We have used two methods to perform 
the fits. First we taken the median value of the emission mea- 
sure in each temperature bin from each Monte Carlo simula- 
tion and fit the resulting distribution. We have also fit each 
distribution individually. Both methods yield consistent re- 
sults. The values for a are indicated on each plot as well as 
in Table 1. The uncertainties indicated in the plots are the 1- 
a standard deviations in the indexes determined from fitting 



The same as Figure 6 but for regions 13-15. Two emission measure distributions are shown for region 13. 
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FIG. 10. — An alternative emission measure distribution derived for region 
1 of the 201 1 April 15. This DEM was computed using a much smaller bin 
size (A log T e = 0.02) than those shown in Figures 6-7 (A log T e = 0.05). 

this would depend on having accurate methods for extrapolat- 
ing the photospheric field into the corona and such extrapola- 
tions have proven difficult to achieve (De Rosa et al. 2009). It 
seems likely that trends observed Figure 9 would also be evi- 
dent in a plot of EM as a function of footpoint field strength, 
but this has yet to be demonstrated. 



FIG. 9. — The total inter-moss emission measure at "warm" temperatures 
(logTe = 6.0-6.2) and hot temperatures (logT e = 6.4-6.8) as a function 
of the total unsigned magnetic flux. The bottom panel shows the power law 
index (a) on the EM distribution between log T e of 6.0 and 6.6. 

each distribution and suggest uncertainties of 10 to 20%. In 
this sample 11 of the 16 EMs have a > 3. However, we also 
measure 5 temperature distributions that are much shallower, 
with a ~ 2, which is similar to the results from Tripathi et al. 
(201 1). It is clear that these shallow EMs are much more com- 
mon in the active regions with the weakest magnetic fields. 

Inspection of the emission measure distributions reveals an 
unexpected trend in the amount of 1 MK emission in the core 
of an active region. In regions 1 through 5 the emission mea- 
sure near 1 MK is often between 10 26 and 10 27 cm~ 5 . In the 
regions with the strongest magnetic fluxes (regions 10-15) the 
emission measure appears to be somewhat smaller, typically 
between 10 25 and 10 26 cm~ 5 . To quantify this we sum the 
emission measure between log T e of 6.0 and 6.2 and plot it as 
a function of total unsigned magnetic flux. As is indicated in 
Figure 9 the emission measure at these lower temperatures is 
inversely proportional to the field strength. This clearly evi- 
dent in Figures 1-3, which show relatively few loops in the 
inter-moss regions in the 171 A for the largest values of mag- 
netic flux. 

The emission measure at the highest temperatures, as ex- 
pected, rises with increasing total unsigned magnetic flux. 
This is also shown in Figure 9. It is important to recognize 
that this comparison between the properties of the inter-moss 
DEM and the total unsigned magnetic flux is not ideal since 
we are comparing an apex property of selected loops with the 
magnetic properties of the entire active region. As pointed 
out by Schrijver (1987), much of the increase in the total un- 
signed magnetic flux simply reflects an increase in the area 
of the active region. The mean field strength also rises with 
increasing active region area, but weakly (also see Fludra & 
Ireland 2008). Ideally we would compare the properties of the 
DEM with the magnetic properties at the loop footpoints, but 



3. DISCUSSION 

We have presented the calculation of emission measure dis- 
tributions for 15 active region observations spanning almost 
an order of magnitude in total unsigned magnetic flux. This 
analysis suggests that the shape of the emission measure dis- 
tribution depends on the magnetic properties of the active re- 
gion. For regions with appreciable magnetic flux the emission 
measure distribution is often strongly peaked at a temperature 
of about 4 MK. For lower levels of magnetic flux, however, we 
do observe shallower temperature distributions. This suggests 
a possible resolution of the varied results presented previously 
(Tripathi et al. 2011; Viall & Klimchuk 2011; Winebarger 
et al. 201 1; Warren et al. 201 1). 

These results are difficult to reconcile with the Parker 
nanoflare model (Parker 1988), at least as it has often been 
interpreted (Cargill 1994; Klimchuk & Cargill 2001; Cargill 
& Klimchuk 2004). As mentioned previously, hydrodynamic 
simulations suggest much flatter emission measure distribu- 
tions than we observe in most of these active regions (Mulu- 
Moore et al. 2011). In the simulations the steepest slopes 
(2.0 < a < 2.3) are obtained for radiative losses based on 
coronal abundances. For all of the inter-moss regions that we 
considered the intensities of the S and Ar emission lines com- 
puted from the DEM are consistent with what is observed, 
indicating that our assumption of coronal abundances is cor- 
rect. It is possible, however, that some of the assumptions 
made in the hydrodynamic simulations, such as constant loop 
cross section or the highly simplified chromosphere, produce 
misleading results. 

It seems likely that high frequency heating that is concen- 
trated at low heights in the solar atmosphere will be able to 
account for the active region properties that we present here. 
The wave heating model described in van Ballegooijen et al. 
(201 1) and Asgari-Targhi & van Ballegooijen (2012) appears 
to be a viable candidate. Detailed simulations, however, are 
required to establish this. 

A number of previous studies have suggested that emission 



A Systematic Survey of High Temperature Emission 



11 



f 












A. 






# 


AIA/I31A 


AIA 4 171 A 


AIA 94 T 


AIA 1 9Tk 


AIA 21 1^ 











> 




AIA 94 A Observed 


AIA 94TEstimated 


/ 

Difference 




10 20 30 40 
AIA 193+171 A Intensity / 116.54 



50 



FIG. 1 1 . — Time-averaged AIA images of a bright point and the quiet Sun. Each image is an average of 300 co-aligned exposures. These data were taken 2010 
March 22 from 12 to 13 UT. The bottom panels show the observed AIA 94 A image compared with the intensities inferred from 171 and 193 A. The differences 
between these images is also shown. The far right panel shows a polynomial fit to the intensities. 



measure analysis is of little utility since the inversion of equa- 
tion 1 is ill-posed (e.g., Craig & Brown 1976; Judge et al. 
1997). It is clear, however, that the general properties of 
active region temperature structure can be determined from 
the available data. We can, for example, safely conclude that 
the emission measure near 4MK is approximately 100 times 
larger than the emission measure near 1 MK in many of these 
active regions. This result is evident in all of the Monte Carlo 
runs and in many different active regions, so is robust against 
perturbations in the observed intensities. It is also clear that 
the detailed structure of the emission measure distributions 
is much more difficult to determine with confidence. Small 
changes in the parameters used in the inversion can lead to 
different results (e.g., Landi et al. 2012). If we run the MCMC 
code with a smaller temperature binning, for example, we ob- 
tain distributions with much more structure. In the example 
shown in Figure 10, the general trend is preserved, but the 



emission measure distribution appears to break up into a se- 
ries of nearly isothermal components (see Landi & Feldman 
2008 for a similar result). Understanding the detailed struc- 
ture of the emission measure distribution will require more 
detailed mathematical analysis. At present, however, develop- 
ing models of the coronal heating process which make predic- 
tions comparable to the observations described here is likely 
to lead to the most rapid progress on this long standing prob- 
lem in solar physics. 

This research was supported by NASA. Hinode is a 
Japanese mission developed and launched by ISAS/JAXA, 
with NAOJ as domestic partner and NASA and STFC (UK) 
as international partners. It is operated by these agencies in 
co-operation with ESA and NSC (Norway). HPW benefited 
greatly from discussions at a International Space Science In- 
stitute meeting on coronal heating lead by Steve Bradshaw 
and Helen Mason. 
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FIG. 12. — An example calculation of the Fe XVIII 93.92 A intensity. The 171 A and 193 A images are combined and scaled to estimate the warm contribution 
to the 94 A channel. The warm emission is then subtracted from the observed 94 A channel image to yield an estimate of the Fe XVIII 93.92 A intensity. The 
bottom panels show the intensity as a function of position across the image. The intensities are averaged along the region indicated by the dotted lines. These 
data are from 201 1 November 8 near 16:34 UT. 

APPENDIX 

AN EMPIRICAL CORRECTION TO THE AIA 94 CHANNEL 

The AIA 94 A channel is contaminated by "warm" emission formed at temperatures much less than the 7.0 MK temperature 
characteristic ofFe XVIII. To illustrate this we have taken one hour of AIA observations (2010 March 22 1 2-1 3 UT) and computed 
time-averaged images from all of the available data. These data were chosen because they contained a large bright point in 
addition to the quiet Sun and show a relatively large range of intensities. The bright point, however, is unlikely to contain any 
Fe XVIII, which would complicate the analysis. The averaging naturally leads to some smearing of the images but is necessary 
to improve the signal to noise. The averaged images for 5 wavelengths are shown in Figure 1 1 . 

Inspection of these images suggests that the warm emission is closest to 193 in morphology. Note the strong contrast between 
the bright point and the quiet Sun, for example. A detailed examination of the loops around the bright point indicates that there 
is also a contribution from cooler emission similar in temperature to 171 A. See (Testa & Reale 2012) for a discussion of stellar 
observations of the this wavelength range. To estimate the intensities in the 94 A channel we consider a polynomial fit to an 
mixture of 171 and 193 A images 

7/m + (l-/)W" 



0.39 

i=0 



116.54 



(Al) 



where the scaling factors derived from the median intensities (1 16.54 and 0.39) have been introduced for convenience. We have 
determined that for / = 0.31 the estimated intensities are closest to what is observed. For this value of / the coefficients to the 
polynomial fit are -7.31 x 10~ 2 , 9.75 x 10~\ 9.90 x 10~ 2 , and -2.84 x 10" 3 . Since there is very little data for very high 
intensities in these data we limit the value of the composite 171/193 A intensity to 30 in using the polynomial fit. The observed 
and estimated 94 A intensities for these data are shown in Figure 1 1 . 

An example set of images is shown in Figure 12. These data were considered by Teriaca et al. (2012) and compared with 
spectroscopic observations of Fe XVIII 974.86 A. Note that this procedure will not work during a flare since Fe XXIV 192.04A 
is likely to contribute to the 193 A channel. This approach will also run into problems for very bright 1 MK emission, such as is 
found in the moss. 

A similar method for isolating the Fe XVIII in the AIA 94 A channel was considered by Reale et al. (2011). They used only 
171 A, however, which does not approximate the contaminating emission as well as a combination of 171 A and 193 A. An 
innovative technique for visualizing the relative contributions of active region emission at various temperatures, including the 
very high temperature Fe XVIII emission, has been presented by Testa & Reale (2012). 
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